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Abstract. We explain the main concepts centered around Sharafutdinov's ray 
transform, its kernel, and the extent to which it can be inverted. It is shown how 
the ray transform emerges naturally in any attempt to reconstruct optical and 
stress tensors within a photoelastic medium from measurements on the state of 
polarization of light beams passing through the strained medium. The problem of 
reconstruction of stress tensors is crucially related to the fact that the ray transform 
has a nontrivial kernel; the latter is described by a theorem for which we provide 
a new proof which is simpler and shorter as in Sharafutdinov's original work, as 
we limit our scope to tensors which are relevant to Photoelasticity. We explain how 
the kernel of the ray transform is related to the decomposition of tensor fields into 
longitudinal and transverse components. The merits of the ray transform as a tool 
for tensor reconstruction are studied by walking through an explicit example of 
reconstructing the eras-component of the stress tensor in a cylindrical photoelastic 
specimen. In order to make the paper self-contained we provide a derivation of 
the basic equations of Integrated Photoelasticity which describe how the presence of 
stress within a photoelastic medium influences the passage of polarized light through 
the material. 
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1. Introduction 

A transparent optical medium which initiahy is isotropic and homoge- 
neous with respect to light propagation may become birefringent when 
subjected to external strain, an effect which is known as Photoelastic- 
ity [1, 2]. The spatially varying tensor of refraction then reflects the 
presence of stresses within the material, opening up the possibility 
of examining internal stresses by means of the change in the state 
of polarization of light beams passing through the strained medium. 
The reconstruction of local optical and stress tensors from the set of 
data collected for all possible directions, and locations, of light beams 
passing through the specimen is called Integrated Photoelasticity [3]. 
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Numerous efforts liave been made in tlie past to tackle tlie problem 
of Integrated Photoclasticity [3-12]. While the reconstruction of the 
optical tensors in the two-dimensional problem is well-established (for 
an overview see [5]), the three-dimensional (3D) case so far has been 
solved only for special cases involving a priori assumptions about the 
symmetry of the 3D stress distribution, the rotation of the principal 
axes of the stress tensor, or the strength of the anisotropy of the di- 
electric tensor. The most general case of 3D tensor tomography with 
arbitrarily shaped bulk matter and no restriction on the degree of 
anisotropy is unsolved. The source of difficulty is the fact that the 
Radon transform which provides a well-established formalism for the 
reconstruction of scalar fields can no longer be utilized if tensor fields 
arc involved. 

In [13], Sharafutdinov has proposed a generalization of the Radon 
transform to the case of symmetric tensor fields of arbitrary degree m 
defined on a, possibly curved, Riemannian manifold of arbitrary di- 
mension n; he has termed his construction the ray transform. Amongst 
many other aspects, his work encompasses the examination of the for- 
mal structure of the ray transform, its relation to the Fourier transform, 
its non-trivial kernel and the problem of inversion of the ray transform 
given the fact that the kernel is non-zero. Sharafutdinov's contribu- 
tion undoubtedly opens up the right path to the goal of solving the 
general problem of Integrated Photoelasticity; however, since his work 
aims to tackle the most general case of symmetric tensors of arbitrary 
degree m in n dimensions, his formalism necessarily exhibits a degree 
of complexity which may be undesirable for practitioners who wish 
to put his theory into a specific physical or engineering context. It 
is therefore an important task to lay out his results in a way which 
focuses on the special case of symmetric {m = 2) tensor fields, and 
vector (m = 1) fields, defined in a three-dimensional Euclidean space 
which is to be identified with the bulk material of the medium, and to 
study the merits of the ray transform as a tool for tensor reconstruction 
within the specific framework of Photoelasticity. 

In the present paper we have attempted to provide some steps in 
this task: We explain the main concepts related to the ray transform, 
its kernel, and the extent to which it can be inverted, for the special 
case of vector- and symmetric tensor fields in IR^. Our goal is on the 
one hand to provide an overview of the mathematical structure of 
the ray transform, and on the other hand to illustrate the method 
by discussing a simple application within the field of Integrated Pho- 
toelasticity. With these objectives in mind we have performed a new 
proof on one of the central theorems in Sharafutdinov's theory, namely 
the kernel of the ray transform; this theorem provides essential insight 
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into the degree to which optical or stress tensors can be reconstructed 
from photoelastic data. Our proof focuses on vector- and symmetric 
2-tensors only, which moreover are reconstructed plane-wise, so that 
they are effectively defined on a two-dimensional R^. This allows us to 
take alternative routes in the various steps of the proof which are not 
available in the most general case. This proof is given along the way 
of explaining the emergence of the ray transform in photoelastic stress 
reconstruction. 

The plan of the paper is as follows: In section 2 we derive the basic 

equations of Integrated Photoelasticity governing the evolution of the 
components of the electric field vector of a polarized light beam passing 
through a medium with spatially varying dielectric tensor. In section 3 
we show how these equations can be used to describe the impact of 
stress in a transparent material on polarized light. In section 4 wc show 
how the concept of the ray transform emerges as we try to reconstruct 
optical and stress tensors from photoelastic measurement data. In sec- 
tion 5 we characterize the kernel of the ray transform of vector and 
symmetric tensor fields, providing a new proof, albeit with a narrower 
scope, to Sharafutdinov's general theorem. In section 6 we show how 
the kernel of the ray transform is related to the decomposition of tensor 
fields into longitudinal and transverse components and explain how this 
knowledge can be utilized to reconstruct stress tensor components in 
a photoelastic specimen. Finally, in section 7 we speculate about the 
extension of the formalism at hand into the high-stress regime using the 
Newton-Kantarovich method. In section 8 we summarize our results. 



2. Basic equations of Integrated Photoelasticity 

We now derive the basic equations of Integrated Photoelasticity un- 
der the straight-line assumption, i.e., light rays passing through pho- 
toelastic media are assumed to propagate along straight lines. This 
assumption is justified as long as anisotropy is sufficiently weak. 

We wish to examine the propagation of a plane polarized electro- 
magnetic wave through a weakly birefringent photoelastic medium. 
The birefringence, i.e., the anisotropy of the dielectric tensor eij,i,j = 
1,2,3, is assumed to be entirely resulting from the stress within the 
material, such that, in the unloaded state, the material is homogeneous 
and isotropic with a real dielectric tensor eij = e 5ij, and e = const. The 
material is assumed to be unmagnetic with vacuum permeability /xq, 
and non-absorbing at least in the optical range of wavelengths, so that 
the material is transparent to visible light. The permittivity of vacuum 
will be denoted as eo- 
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We start with the source-free Maxwell equations in the dielectric 
medium, 

VxH = D , VxE = -B , (la) 

V«E = , V«B = , (lb) 

where we have ignored the microscopic sources of the dielectric; they 
are taken into account phenomenologically by the spatially varying di- 
electric tensor e^j , which links the electric field E and the displacement 
field D according to 

D, = e^jEj , (2) 

where the convention of summing over double indices is adopted. Fur- 
thermore, the effective charge density due to polarization effects is 
ignored, 

Peff = -V«P = . (3) 

If the second equation in (la) is inserted into the first, we obtain the 
wave equation 

-AE + /xo D = , (4) 

where A is the three-dimensional Laplace operator in Cartesian coor- 
dinates. 

We now assume a harmonic time dependence ~ e^*'^* of all fields, 
and a wave propagation in the z-direction. We may neglect the x- 
and y-dependence of the fields, since the electromagnetic energy in the 
geometrical-optics limit propagates along narrow tubes of light rays in 
such a way that, along the cross-section of a given tube, the fields may 
be treated as independent of x and y. Then (4) becomes a Helmholtz 
equation 

^+A^-'D = . (5) 

Provided that birefringence is weak, i.e., the load on the material is not 
too strong, the component of the electric field in the direction of the 
wave vector can be neglected. This leads to an Ansatz for the electric 
field as follows: 

E(z,i) = ^He{^i(z) 61+^2(^)62} e^'^"-^"^* , (6a) 
k = — , u = . (6b) 

This Ansatz contains the assumption that the phase velocity u of the 
wave is homogeneous and isotropic, and is determined by the permittiv- 
ity e of the unstressed medium; this is true under the same conditions 
under which the z-component of the electric field can be neglected, i.e. 
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sufficiently weak anisotropy. Eqs. (6) express the straight-line assump- 
tion mentioned at the beginning of tliis section. 

Now the the second derivative in z of the complex electric field in 
(6) is 

E'/ = A'^e'''' + 2ikA'ie'''' -k'^Aie'''' , i = l,2 , (7) 

where a prime indicates a derivative d/dz. If the medium were ho- 
mogeneous, the amplitudes were independent of z, just as within 
vacuum. The inhomogeneity of the medium makes them ^-dependent, 
but as long as birefringence is weak, this dependence is sufficiently weak 
to neglect the second derivative A" in (7). If we then insert (7) into (5) 
and divide by 2ik, we obtain 

We introduce a constant 

C„ = |L = ;^ , (9) 

with the help of which (8) can be written as 
dA- 

-^=A'i=iCo {eij Aj-eAi) , i = l,2 , (10a) 
or as a matrix equation 

(10b) 

The complex two-component vector A = (^41,^2)^ is called the Jones 
vector of the light beam. Eq. (10b) can be converted into an operator 
equation as follows: We may partition the photoelastic medium into 
many thin slices perpendicular to the direction of the wave propagation 
e^. Each slice acts on the Jones vector of the light beam either as 
a retarder, i.e., introducing a phase shift between two orthogonally 
polarized components of the wave, or as a rotator, i.e., rotating the 
plane of polarization of a plane-polarized wave, or both. This implies 
that each slice acts on the local Jones vector as a unitary matrix, and 
hence the same must be true for the total bulk of the photoelastic 
specimen. If Zi is some point before the medium (possibly Zi = —00) 
and Zf is some point behind the medium (possibly Zf = -|-oo), then the 
Jones vectors Ai and Af at these points must be related by a unitary 
transformation U{zf, zi). In fact, such a transformation must exist for 
any pair of points z, z' , so that we can write 

A^ = U{z, z') A^, for all z' . (11) 
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The matrix U{z, z') must satisfy an equation analogous to (10b), 

■^U{z,z') = iM{z)U{z,z') . (12) 
The matrices U{z, z') obviously must have the group property, 

U{z3,Zi) = U{Z3,Z2) U{z2,Zi) , 

U{ZI,Z2) = U{Z2,ZI)-' , U{Z,Z) = l2 . ^ ' 



3. Stress-optical relations 

The stress- optical equations linking the dielectric tensor e^- with the 
stress tensor aij, 

eij = e6ij + Ciaij + C2ti(Tdij , (14) 

were discovered long ago by Maxwell [14] . Since the material is assumed 
homogeneous with respect to its stress-optical properties, the quantities 
Ci and C2 are constants, characterising the stress-optical properties of 
the specimen. In (14), Sij and aij depend on the spatial coordinates 
{x,y,z), but e is constant throughout the material, as mentioned in 
section 2. Using (14) we can express the equations (10) in terms of the 
stress tensor aij rather than the dielectric tensor eij. We start again 
with (10b) and substitute (14) for the components of the dielectric 
tensor in the matrix M{z); this gives 

where the trace tr a runs over all three indices, 

3 

tTa=Y^ann ■ (16) 

n=l 

We now perform the split 

Ci (Til + C2 tr (7 = ^ Ci {an - (T22) + B , (17a) 
Ci £722 + C2 tr (7 = Ci (an - (722) + B , (17b) 

where 

B= - Ci{au+CT22) + C2tTa . (17c) 
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When inserted into (15), the matrix M can be written as 



M{z) = CoCi 



i((Jii-C722) 0"i2 



2 
0'21 



-i((7ll - CJ22J 



N{z) + C0BI2 



+ Co S I2 = (18a) 
(18b) 



We sec that the term CqB just introduces a phase common to both 
components Ai and A2 of the Jones vector. This phase is physically 
immaterial and can be absorbed by defining a "reduced" Jones vector 
Aj., 



Ar = A exp I -iCo J dz B{z) 



(19) 



The dynamical equations (10) expressed in terms of Aj. then take the 
simpler form 



^^Mz) = iN{z)Ar{z) , 



N{z) 



— Co Ci 
2 



Cll — (722 2(712 

2(721 — (0"11 — 0-22). 



(20a) 



(20b) 



Since Jones vectors at different points are unitarily related, see (11), 
eq. (20a) gives rise to an operator differential equation analogous to 
(10b), 

^U{z,z') = iN{z)U{z,z') , U{z,z) = 12 . (21) 
dz 

This can be converted into an integral equation 



U{z, z') = l2 + i j dzi iV(zi) U{zi, z') . (22) 

z' 

A formal solution is given by the Bom-Neumann series 

z z z\ 

U{z,z') = I2 + y dzi iN{zi) + j dzi iN{zi) j dz2 iN{z2) + ■■■ . 

z' z' z' 

(23) 

If the specimen is located such that the ^;-axis intersects the medium, 
and z, z' refer to points before and behind the medium, respectively, 
then the matrix U{z,z') contains all the data acquired from mea- 
surements of the phase retardation between orthogonal components of 
plane-polarized light, and the rotation of the plane of polarized light, 
along that line of intersection. 
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For weak birefringence, (23) can be truncated after the first-order 
term, 

z 

U{z,z') = l2 + i j dziN{zi) . (24) 

z' 

Since the medium occupies a bounded region in space, the hniits z and 
z' in the integral in (24) may be taken as ±00. On inserting (20b) into 
(24) we see that the evolution of the state of polarization along the 
light ray is determined by the two integrals 

00 

Iw = J dz (T12 , (25a) 

—00 
00 

lu = j dz ((Til — 0-22) ■ (25b) 

—00 

The significance of the notation Iw,Iu will become clear shortly. 



4. The ray transform in Photoelasticity 

We now wish to examine the stress within a cylindrical object which 
occupies a cylindrical domain G = D x (a, 6), where D is a two- 
dimensional region in the xy-p\ane and the boundary dD of D is a 
strictly convex smooth curve. By B = dD x (a, b) we denote the lateral 
surface of the cylinder G. The stress tensor is supposed to be smooth 
in G and on B, and satisfies the equilibrium conditions 

everywhere. Furthermore, the absence of external forces on the lateral 
boundary dB implies that 

aij Hj = , (27) 

where (ni, n2, 0) is the unit outward normal to dB on the lateral surface 
B of the cylinder. Strain is therefore applied only at the top (z = b) 
and the bottom {z = a) of the object. 

It is now assumed that measurements analogous to (24) are per- 
formed along all horizontal straight lines in IR^ defined hy t ^ X + 
t^{a), where X = (x,y,z), {x,y) varies through the points in the xy- 
plane, z takes values in the interval (a, b), the angle a varies in (0, 27r), 
and 

^(a) = (cosa, sina, 0) (28) 
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is the horizontal unit vector in the direction of the hne. These hnes 
generalize the measurements along the z-axis which were discussed in 
section 2. For the straight line passing through the point X along the 
direction ^ we now choose an adapted right-handed coordinate system 
(r), z, t) in which the direction of propagation of the light beam is along 
the positive t-axis (the 77-coordinate line is perpendicular to the t- 
coordinate line, and the new z-direction coincides with the old one). 
On performing a measurement along the line so defined we then obtain 
quantities analogous to (25), 

00 

Iw{X,^) = J dtar,z , (29a) 

—00 
00 

Iu{X,0 = J dt {arjr, - a,,) , (29b) 
—00 

the notation on the left-hand side indicating that these lines pass 
through the point X and are directed along the unit vector ^. The 
tensor components in these integrals can be expressed in terms of the 
laboratory-frame components aij as 

(o-jjj? - (Tzz) {t, 0, 0) = sin^ a {an - (733) (^ + ^0 - 
— 2 sin a cos a a\2 i?^ + ^0 + 

-h cos^ a (0-22 - 0-33) (X tCi , (30a) 
(7^^(t,0,0) = -sinac7i3(X-Fte) + 

cos a 0-23 (X + ^0 > (30b) 
(7zz(t,0,0) = cT33(X + tO • (30c) 

We now define a vector field u; in IR^, where 

= I (x, y, z) G (x, y) G , z e (a, 6) fixed} , (31) 

by 

{w\,W2) = {o^-i.-on) , (32a) 
and a symmetric tensor field u by 

(nil, U12, ^21, U22) = (0-22 - o"33,-o'i2,-cr2i,o-ii - o"33) • (32b) 

If (30) is expressed in terms of the components of the tensors u and w 
and inserted into (29) we obtain 

Iw{X,^) = j dtwm^m , (33a) 

Iu{X,0 = Jdt U^n ■ (33b) 
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According to [13], the collection of quantities Iw{X,^), taken for all 
X € IR^ and ^ = ^(a) with a € [0, 27r], constitutes the ray transform of 
the two-dimensional vector field w; whilst the collection of all Iu{X, S^), 
with the same range of variables, constitutes the ray transform for the 
two-dimensional symmetric tensor u. 

The general definition of the ray transform of a symmetric tensor 
field Tij^...i,n of degree m defined in R"-, which includes both Iw and lu 
in (33) as special cases, is given by 

ITiX, 0= jdt Ti„„i^ {X + to • ■ ■ ^im , (34) 

where X,^ G IR*^ with |^| = 1, i.e. ^ lies on the unit sphere Sn-i in 
IR". Since the ray transform is constant for all X' lying on the line 
t^X + t^, 

ITiX + T^,0=IT{X,0 , (35) 
we can restrict the range of the variable X to the subspace 

C± = {x G R"| } (36) 

of IR" orthogonal to ^, since any component of X in the direction of ^ 
may be set to zero. It follows then that the ray transform is really a 
function on the tangent bundle TSn-i of the unit sphere Sn-i, 

r5„_i = {(x,0eR"xS„_i|x±O ■ (37) 

As a consequence, the Fourier transform with respect to X may be 
restricted to the subspace 

^^^'^^ = (27r)(»-^)/^ / ^^""'(-^^) J^(x^Oe-'<^■^'> , kec± . 

(38) 

The Fourier transform (38) of the ray transform is related to the 
Fourier transform T of the tensor T as follows: 

/r(fc,0 = (27r)V2f,,...^(fc)Cn...Ci. , ke^i_ . (39) 
This can be proven easily: Inserting (34) into (38) we find 

^(^^ = ^^Jn-iy2 J dV--\X') J dt r,...,^ • ■ ■ (ir. e-^^'^'^'^ . 

(40) 

Since |^| = 1 and X'_L^, we can introduce coordinates X on IR"' by 
X = {X',t), in which case dV-^{X')dt = dV{X). Furthermore, 
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since k-L^ by definition (38) of the Fourier transform, the argument of 
the exponential is 

{k,X) = {k,X' + tO = {k,X') . (41) 
If these relations are inserted into (40), (39) follows. 

5. The kernel of the ray transform 

Our task is now to reconstruct the components of the stress tensor from 
the ray transforms (33) as far as possible. We shall sec that, for the 
given experimental setup, involving horizontal light rays only, we can 
only retrieve one component of the stress tensor, namely (733. In order to 
proceed we need to understand that the ray transform, regarded as an 
operator / sending symmetric tensor fields T (including fields of degree 
m = 1, i.e. vector fields) T to the set of quantities IT, has a nontrivial 
kernel: Given two tensor fields T,^}^ and T^^^ , does /T^^) = /T^^) 
imply T(^) = T(2)? The answer is no, as proved in a theorem in [13]. 
The original version of this theorem is designed to encompass the most 
general case of symmetric m-tensors defined on IR", with no reference 
to a specific physical application or context. Accordingly, the proof, 
encumbered with this degree of generality, occupies a substantial num- 
ber of pages. However, this generality is not required for the particular 
purposes of reconstructing vector- and symmetric 2-tensor fields within 
the context of Photoelasticity. If we limit our scope to these types of 
fields, the proof of the theorem can be performed along alternative 
routes which are shorter than in the original work. In this section we 
provide this new proof for the two special cases of a vector field w and 
a symmetric tensor field u defined in IR^ as discussed in (32, 33): 

THEOREM 5.1 (Kernel of ray transform). 

A. Let w be a smooth vector field in IR^ with compact support. Then 
the following statements are equivalent: 

a) Iw = . 

b) There exists a compactly supported scalar field (p on IR^ such 

that its support is contained in the convex hull of the support 
ofw, and 

Wi = di4> , i = l,2 . (42a) 

c) The identity 

diW2 - d2Wi = (42b) 

holds. 
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B. Letu be a symmetric tensor field on with compact support. Then 
the following statements are equivalent: 

a) lu = . 

b) There exists a compactly supported vector field v in R? such 

that its support is contained in the convex hull of the support 
of u, and 

'^ij = \ + ' « = 1, 2 . (42c) 

c) The identity 

dfu22+diuu-2did2Ui2 = , (42d) 

holds. 



Proof. 
5.1.1. Case A. 

The equivalence (A.b) <^ (A.c) is well-known in standard Vector Anal- 
ysis, where (42b) expresses the fact that the vector field Wi must be 
curl-free in order for a potential (f) to exist. We therefore shall not prove 
this equivalence here. 

Wc now prove (A.b) =^ (A.a): Inserting (42a) into the first integral 
of (33) yields 

IwiX, 0= j dt Udm(f>{X + tO = Jdtj^(l>{X + tO = , (43) 

since (p has compact support and therefore must vanish at t = ±00. 

Finally we prove (A.a) =^ (A.c): To this end we observe that the 
function diW2 — (?2tfi is a scalar with respect to rotations in the plane 
R^. This is obvious if we extend the field Wi,i = 1,2, smoothly to a 
three-dimensional vector field = 1,2,3, defined in R"^, such that 
wi = wi and W2 = W2 on the plane z = const. Then (42b) is just 
the z-component of the curl of w, which transforms like a scalar under 
S'0(2)-rotations in this plane. It follows that the collection of integrals 

j dt {diW2 - d2Wi)iX + to , (44) 

taken for all directions ^ in IR^, is just the two-dimensional Radon 
transform of this scalar field. If all of the integrals (44) were zero, the 
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invertibility of the Radon transform would imply that the integrand 
must vanish and (42b) must hold. But this is indeed the case: Let us 
take the derivative of the equation Iw = with respect to xi, 

d\Iw = J dt {cos 6 wi^i + sin 9 W2,i} = , (45) 

where we use the abbreviation wi^i = diwi, etc. On account of the fact 
that wi has compact support and therefore must vanish at infinity we 
have ^ 

J dt {cos9wi^i + sin. 9 wi^2} = J ~ ^ ' 

so that (45) gives 

dilw = sine j dt{w2,i - ^1,2} = . (47) 

From d2lw = we infer a similar equation, with sinO replaced by cos 9 
in front of the integral. These equations must hold for all 9, so that the 
line integrals (44) all vanish. This proves the statement. 

This finishes the proof of case (A.). 
5.1.2. Case B. 

Our strategy here is different, since we cannot employ standard tech- 
niques from Vector Analysis. 

We first prove the equivalence (B.b) (B.c): The implication (B.b) 
(B.c) is trivial, and follows immediately by inserting (42c) into (42d). 

Now we prove (B.c) =^> (B.b): To this end we define a two-dimensional 
vector field {vi,V2) as follows: 

2-K 

Vi = ^jd9 j dtUmn{X + ti)Uin^i , (48) 
-00 

where i = 1,2, ^1 = cos^, ^2 = sin^. We can now compute the 
expression didivi = by using the fact that 

^ uii = cos ^ uii^i -I- sin ^ nii,2 • (49) 
Making use of the assumption (42d) we are able to derive 

■"1,11 = . (50a) 

In a similar way we find 

-"1,12 = wii,2 . (50b) 
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Eqs. (50) imply that divi = uu + c, where c is a constant. However, 
uii has compact support and therefore vanishes at infinity, whilst vi 
by construction behaves like 1/|X| for \X\ oo and, in particular, 
vanishes at infinity . It follows that c = 0, and hence 

Uu = diVi , (51) 

which proves eq. (42c) for the case i = j = 1. The proof for Ui2 and 
U22 proceeds in exactly analogous a manner. 

Next we prove (B.b) =^ (B.a): Inserting (42c) into (33b) yields 

luiX, = fdt UindmVn{X + 1 = J ^Hn VniX + tO = , 

(52) 

since v has compact support. 

Finally we prove (B.c) (B.a): Using the definition (33b) of lu, we 
derive the equations 

(/u),ii = j dtly cos^6''Uii,ii + 

+ 2sin0cos^'Ui2,ii + sin^ 0U22,ii| , (53a) 

{lu), 12 = j dt{ COS^6''Uii,i2 + 

+ 2sin6'cos6''Ui2,i2 + sin^ 6'u22,i2| , (53b) 

(7u),22 = j dt { COS^ Uu,22 + 

+ 2 sin 6* cos 9 ^12,22 + sin^ 6 ^22,22 1 • (53c) 

Starting with (53b) we find on using (42d) that 

(-f ii),i2 = j dt^cosB^ iiii,2 + sin 6* ^ ^22,1 | = . (54) 

Using the same technique we obtain for (53a, 53c) 

(/n),n = (/n),22 = . (55) 

It follows that 

(/U),i = /i(^) , (/u),2 = /2(^) . (56) 

The functions /i and /2 can be evaluated at points X arbitrarily close 
to infinity, so using the fact that u has compact support implies that 
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/i = /2 = 0. Then, repeating the same argument for lu^i and Iu^2 
shows that we must have lu = 0. 

This finishes the proof of case (B.). 

This finishes the proof of theorem 5.1. ■ 

From eqs. (42b) and (42d) we learn that elements of the kernel of 
the ray transform can be characterized by the identical vanishing of 
a differential expression of the vector /tensor components Wi,Uij. The 
differential operator acting in (42b) and (42d) is generically called the 
Saint- Venant operator [13]. Its vanishing can be understood as an inte- 
grability condition ensuring the existence of the scalar/vector fields (f> 
and V occuring in cqs. (42a) and (42c). Wc sec that in the simplest case 
of a vector field w, the condition for the vanishing of the ray transform 
is equivalent to the statement that w is the gradient of a potential 
(f), and (42b) is nothing but the necessary and sufficient condition for 
the existence of the potential. The case of a symmetric m-tensor in 
provides a nontrivial generalization of this scenario. 



6. Longitudinal and transverse tensor components 

The significance of the kernel of the ray transform is further eluci- 
dated by considering the standard decomposition of symmetric tensor 
fields of degree m (including vector fields) into transverse Tj_i^ i^ and 
longitudinal {dv)i^,„i^ components, 

Ti^...im = r±ii..i™ + {dv\...i^ . (57) 

The transverse part rj_i^...i^ is characterized by its vanishing diver- 
gence, di^ Tj^n^.i^ = 0, while the longitudinal part dv is the sym- 
metrized covariant derivative of the symmetric (m — 1) tensor v, and 
is given by 

^'^'")n-im^ ~\ X] ^V(l)-V(m-l),»7r(m) ' (^^) 
III. ^ 

where Sm is the group of permutations of m elements. The transverse 
part of T is determined by the projection of the Fourier transform 
Ti-^...im{^) of ^ii- -«m(^) o^to subspace orthogonal to the direction 
of propagation k of the Fourier mode e*^'^'^^ 

Tu....^{X) = J a!"fcfx,„..,„(fc)e^<'='^> , (59a) 

T±n...ijk) = Xi,j, ■ ■ ■ Xi^j^ fj,.„jjk) , (59b) 
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where 

= % - -gr 

is the projector onto the subspace perpendicular to k. By construction, 
this projector leaves the Fourier transform 2\ of the transverse com- 
ponent Tl invariant, while it annihilates the Fourier transform T — T± 
of the longitudinal component T — T±. 

For our vector field w and tensor field u, the decompositions (57) 
take the special form 

Wi = w±i + di(l) , di w±i = , (60a) 

Uij = u±ij + ^ {di Vj + dj Vi) , dj u±ij = . (60b) 

On comparing (42a, 42c) with (60a, 60b) we see that the kernel of the 
ray transforms Iw, lu consists precisely of the longitudinal components 
d^, dv of the vector/tensor field w,u\ It follows that these components 
cannot be retrieved from the ray transforms of these fields. This result 
is obviously of central import for any attempt to reconstruct optical 
and stress tensors from photoelastic measurements. 

Endowed with this knowledge we can now turn to the question 
to which extent the ray transform may be inverted: To this end wc 
introduce the integral moments {ii"^IT)j^,,,j^{X) of the ray transform 
{IT){X,^) of the symmetric tensor ri^,,,j^(X) of degree m on R", 

{l^'^mh-jm (^) = 77^ / d^n-iii) ■ ■ ■ iir. (IT) (X, , (61) 

Sn-l 

where Sn-i is the unit sphere in R"^ with (n — l)-dimcnsional volume 
con-i and volume element dO„_i(^). For our vector field w and sym- 
metric tensor field u in IR^, the integral moments of the ray transform 
Iw{X,^) are 

{ti'lw)i{X) = ^Jde^i {Iw){X,0 , (62a) 



27r 

[ii^Iu),,{X) = dei,ij{Iu){X,i) . (62b) 



On taking (two-dimensional) Fourier transforms of eqs. (62a) and (62b) 
we arrive after a lengthy calculation at 

{ji^Iw^Xk) = jrr Wm{k) \im{k) , (63a) 
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(li^Iu],Xk) = - Umn{k) \k\~'^ X 
\ y 13 6 



(63b) 



The presence of the projector Ajj in these formulas cancels out the 
longitudinal components of Wm and Umm while leaving their transverse 
components wj_yn and u±rnn invariant. It follows that we can replace 
Wm and Umn by w±m and u±mn in formulas (63). On taking the trace 
of the resulting equations we can then invert them for iv± and u±, 



u±ij{k) = ^ \k\ (/hi),, -\\k\Y. (?^)„„ ■ (64b) 

Eqs. (64) provide the explicit answer to the question to which degree 

the components of a vector field w and a symmetric tensor field u on 
can be reconstructed from their ray transforms Iw and lu: It is pre- 
cisely the transverse components of these fields which can be retrieved 
in full, whereas the longitudinal components, these being symmetric 
covariant derivatives of lower-rank tensors, must remain undetermined. 
Performing the inverse Fourier transform on (64) we find that 

wUX) = [liHw] {X) = 



(64a) 



= U 



H^Iu {X) 



- - 1*1 (J^Tu) . (65b) 

Here, Wj [/U-'^/it;] (X) indicates a linear functional with respect to fi^Iw, 
being parametrized by the spatial coordinates X; and similarly for U. 
The right-hand sides of (65) contain the ray transforms Iw and lu of 
the vector- and tensor field, and can therefore be regarded as known 
functions of the data collected from measurements. Now, from theorem 
5.1 we know that the Saint- Vernant operator (42b, 42d) annihilates 
precisely the longitudinal component of a vector- or symmetric tensor 
field [and therefore must be related to products of the projectors (59c) 
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by Fourier transformation]. Thus, if the Saint- Venant operator acts on 
w and u in the manner of (42b) and (42d), only the transverse compo- 
nents w_i_ and n_|_ of w and u survive. For these transverse components 
we now substitute expressions (65); thus, we obtain quantities 

W = 5i>V2-a2Wi , (66a) 
V = df U22 + diUn-2 di 82 U12 , (66b) 

where W and U represent known functions of measurement data. On 
the other hand, we can insert the values (32a) and (32b) for the com- 
ponents of the vector w and tensor u into the left-hand sides of (65), 
so that 

W = -diais- 82 (723 , (67a) 
V = 8f ((Til - (T33) -I- di ((T22 - (^33) + 2 ^1^2 (712 ■ (67b) 

We now regard (67) as a system of equations for the unknown aij , where 
the left-hand sides W and U arc known through (66); furthermore, the 
equilibrium conditions (26) and the boundary conditions (27) are added 
to this system. The third equation in (26) implies that 

y = 8s (733 • (68) 

Furthermore, differentiating the first equation in (26) with respect to 
xi, the second with respect to X2 and adding the results gives, on using 
(67a), 

SaW-U = (af + ai) (733 ■ (69) 

This last equation defines a Dirichlet problem in IR^ for the unknown 
(733, the left-hand side containing a known function of measurement 
data. 

The value of (733 on the lateral surface B of the cylindrical region G 
is determined by the ray transform lu: Let Iu{X, X') denote the ray 
transform of u along a line which connects two points X = {x, y, z) 
and X' = {x',y',z) at the same height z on the cylindrical specimen, 
where both X and X' belong to the boundary B. By taking the limit 
X' ^ X, the integrand in lu is nonvanishing only on an infinitcsimally 
short interval in the neighbourhood of the point X, which becomes 
tangential to B in the limit. Using the stress-free condition (27) of the 
lateral surface then yields 

at every point X E B. Relations (69, 70) are sufficient to uniquely 
determine the component (733 from the measured data, i.e., W and U, 
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since the value of C733 on the boundary is given by (70), hence the 
Dirichlet problem (69) can be solved uniquely for any given z € (a, 6). 
Alternatively, the Dirichlet problem may be solved for just one € 
(a, 6) and then (68) may be integrated to obtain C733 everywhere. 

Finally, we might wonder whether there is any more information to 
be retrieved from the ray transforms Iw and /u, or alternatively W and 
m, as so far we could only reconstruct one tensor component from these 
data. The answer is no: As proven in [13], if any two stress tensors 

and cr^^^ satisfy (26, 27) and 0-33^ = alg'* then they produce the same 
ray transforms, IvP-^ = Iw^"^^ and lu^^' = Iu^'^\ As a consequence, the 
retrieval of one tensor component is the maximum information to be 
gleaned from a ray transform which is performed for horizontal lines 
only. If we wish to reconstruct other tensor components, the orientation 
of the incident light beam must be changed accordingly. 

7. Newton-Kantarovich method 

Sharafutdinov's theory shows that in the linearised problem of pho- 
toelastic tomography the transverse part of the stress tensor can be 
recovered from complete photoelastic data, and that this inversion 
is stable under suitable smoothness assumptions on a. The forward 
problem represented by the solution of the ODE (21) along each ray is 
non- linear in a, its Frechet derivative at cr = being the ray transform 
(24). Unlike the linearization about non-zero stress, this transform has 
an explicit Fourier inversion formula. It might therefore be possible to 
extend the reconstruction method outlined above into the high-stress 
regime by using the Newton-Kantarovich method with fixed derivative, 
which means that while the forward problem is solved successively for 
each updated stress, the difference between predicted and measured 
optical data is used as boundary data for the same ray transform at 
each iteration. While the convergence of this method is linear compared 
with the quadratic convergence when the derivative is updated, it is 
expected that the efficient computation of the linear step will make the 
former more efficient numerically. These ideas will be investigated in 
upcoming papers. 

8. Summary 

We have described the main ideas centered around Sharafutdinov's ray 
transform which represents a generalization of the Radon transform 
to tensor-field integrands. The tensorial character of the field to be 
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reconstructed manifests in a non-trivial kernel of the ray transform, 
as a consequence of which the ray transform can be inverted only 
with respect to the transverse components of the tensors, while the 
longitudinal components are lost. We have provided an alternative 
proof to the main theorem about the kernel of the ray transform as 
given in Sharafutdinov's original work, which is simpler since the scope 
of our proof is limited to those cases which are relevant for the field 
of Integrated Photoelasticity. In order to illustrate the merits of the 
ray transform as a tool for tensor reconstruction, a simple example 
in Photoelasticity, namely the reconstruction of the (J33 component 
of the stress tensor inside a birefringent medium from its horizontal 
ray transform, has been discussed. The inversion of the ray transform 
for the transverse components of the stress tensor together with the 
equilibrium conditions on aij and the strcss-frcc condition of the lateral 
boundary, resulted in a well-posed Dirichlet problem for the compo- 
nent (T33 of the stress tensor which admits a unique solution. Finally 
we suggested a non-linear inversion algorithm based on the Ncwton- 
Kantarovich method with fixed derivative which might be capable of 
utilizing the ray transform in the high-stress regime. 
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